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The microscopic origin of metamagnetism and metamag- 
netic transitions in strongly anisotropic antiferromagnets is 
investigated within a quantum mechanical theory of corre- 
lated electrons. To this end the Hubbard model with stag- 
gered magnetization m s t along an easy axis e in a magnetic 
field H || e is studied both analytically and numerically within 
the dynamical mean field theory (DMFT). At intermediate 
couplings the self-consistent DMFT equations, which become 
exact in the limit of large coordination number, are solved by 
finite temperature Quantum Monte Carlo techniques. The 
temperature and magnetic field dependence of the homoge- 
neous and staggered magnetization are calculated and the 
magnetic phase diagram is constructed. At half filling the 
metamagnetic transitions are found to change from first order 
at low temperatures to second order near the Neel tempera- 
ture, implying the existence of a multicritical point. Doping 
with holes or electrons has a strong effect: the system be- 
comes metallic, the electronic compressibility increases and 
the critical temperatures and fields decrease. These results 
are related to known properties of insulating metamagnets 
such as FeBr2, metallic metamagnets such as UPdGe, and 
the giant and colossal magnetoresistance found in a number 
of magnetic bulk systems. 

PACS: 71.27+a, 75.10.Lp, 75.30Kz 



I. INTRODUCTION 

While ferromagnetism is known since antiquity, anti- 
ferromagnetism was only discovered in this century. It is 
not widely known that the concept of antiferromagnetic 
order wasproposed independently by Neel [Q in 1932 and 
Landau M in 1933. Both sought to explain the, at that 
time, puzzling low temperature behavior of the magnetic 
susceptibility of certain materials: of metals such as Cr 
and Mn in the case of Neel, and of insulators with layered 
structure such as the chlorides of Cr, bivalent Fe, Co and 
Ni in the case of Landau. While Neel correctly suggested 
the existence of interpenetrating sublattices in Cr and 
Mn with opposite magnetization |IJ , Landau equally cor- 
rectly predicted the existence of stacks of ferromagncti- 
cally ordered layers whose magnetization alternates from 
layer to layer g|. In both cases the total spontaneous 
magnetization adds up to zero. Assuming the interlayer 
coupling to be weak, Landau argued that a relatively 
small magnetic field would be sufficient to modify the 
mutual orientation of the moments in each layer. This 
leads to deviations from the linear dependence of the to- 
tal moment on the field, i.e. to an anomalous increase 



of the susceptibility, and finally - at high fields - to a 
saturation of the magnetization. Such a behavior was in- 
deed observed by Becquerel and van den Handel in 1939 
H in the mineral mesitite (carbonate of Fe and Mg) at 
low temperatures. Not being aware of Landau's or Neel's 
work they could not explain their observation in terms of 
ferro- and paramagnetism, and therefore suggested for it 
the name metamagnetism |4|. 

Qualitatively similar, but even more drastic magne- 
tization effects were later observed in many other sys- 
tems of which FeCl2 and Dy3AlsOi2 (DAG) are well- 
studied prototypes j^j. These materials are insulators 
where the valence electrons are localized at the Fe and 
Dy ions, respectively. The resulting local moments order 
antiferromagnetically and are constrained to lie along an 
easy axis e, implying a strong anisotropy such that a 
spin-flop transition cannot occur. Under the influence 
of a large magnetic field H || e the staggered magne- 
tization vanishes in a first or second order phase tran- 
sition, the so-called metamagnetic transition. Apart 
from the thoroughly investigated materials mentioned 
above, there are also conducting systems in that class, 
e.g. Uranium-based mixed systems H|), SmMn 2 Ge 2 0, 
and TbRlWr^Sia §. 

Today the term "metamagnetic transition" is used in 
a much wider sense PJlOl , namely whenever the homoge- 
neous susceptibility x(H) has a maximum at some value 
H c , with m(H) being strongly enhanced for H > H c . 
Metamagnetism is then found to be a rather common 
phenomenon which occurs also in spin-flop antiferromag- 
nets (e.g. the parent compound of high-T c superconduc- 
tivity, La2Cu04 ]To[|), strongly exchange-enhanced para- 
magnets (e.g. TiBe2, YC02 P[), heavy fermion and inter- 
mediate valence systems (e.g. CeRu2Si2, UPta p0[). 

In this paper we will be concerned only with metamag- 
netism in strongly anisotropic antiferromagnets, several 
of which are known to have a very interesting H—T phase 
diagram {H: internal magnetic field, T: temperature). 
In particular, in the insulating systems one often finds 
a tricritical point at which the first order phase transi- 
tion becomes second order. Theoretical investigations of 
tricritical points began with the work of Landau, who de- 
scribed multicritical behavior within his phenomenolog- 
ical theory of phase transitions Jll|]. Clearly, a genuine 
understanding of the origin of tricritical points etc. re- 
quires investigations on a more microscopic level. In the 
case of strongly anisotropic metamagnets those investi- 
gations where sofar restricted to the insulating systems, 
such as FeBr2 . They are usually based on the Ising model 
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with more than one interaction, in a magnetic field, on a 
simple cubic lattice 0|, e.g. 



H 
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NNN 



where the summations extend over the nearest neighbors 
(NN) and next nearest neighbors (NNN) of every site. 
For J, J' > one has an antiferromagnetic (AF) cou- 
pling between the Z nearest neighbor spins and a ferro- 
magnetic coupling between the Z' next nearest neighbors 
p3| . It was pointed out by Kincaid and Cohen [Q that 
in mean-field theory a tricritical point |l5| ] as in Fig. la 
exists only for R = Z'J'/(ZJ) > 3/5, while for R < 3/5 
this point separates into a critical endpoint (CE) and a 
bicritical endpoint (BCE) (see Fig. lb). The latter be- 
havior, especially the finite angle between the two tran- 
sition lines at CE and the pronounced maximum at the 
second order line, is qualitatively very similar to that 
observed in FeBr2 p6| . However, the first order line be- 
tween CE and BCE has so far not been observed - nei- 
ther experimentally, nor even theoretically when evaluat- 
ing (|l|) beyond mean-field theory p^-pj[| . Most recently, 
by measurement of the excess magnetization and anoma- 
lous susceptibility loss J|(| , the specific heat [^l[ , and by 
neutron scattering (22) a strip-shaped regime of strong 
non-critical fluctuations in the H a — T plane (where H a 
is the applied field) of FeBr2 was reported which is at 
least reminiscent of the critical line CE <-> BCE. This 
unusual behavior was then also found theoretically by 
Selke [jl9| who evaluated (]l|) and also a more detailed 
model [p3[ for weak ferromagnetic coupling (R = 0.4) 
and large Z using Monte Carlo techniques. There does 
not yet exist a microscopic theory for the conducting sys- 
tems such as the Uranium-based mixed systems f|U, 
because this requires a fully quantum mechanical treat- 
ment of itinerant, correlated electrons. First steps in this 
direction, where genuine correlation effects were, how- 
ever, neglected, are the se mi-p henomenological theories 
of Wohlfarth and Rhodes |2j| for metamagnetic phase 
transitions in paramagnets, and of Moriya and Usami 
|p5| for the coexistence of ferro- and antiferromagnetism 
in itinerant electron systems. It is the purpose of this pa- 
per to examine the origin of metamagnetism in strongly 
anisotropic antiferromagnets from a microscopic, quan- 
tum mechanical point of view. To this end we will study 
the Hubbard model in the presence of a strong anisotropy 
direction, with a magnetic field along this direction, since 
it is the simplest microscopic model that describes insu- 
lating and metallic, spin-localized and bandlikc antifer- 
romagnets with easy axis in an external magnetic field. 

In Sec. H the underlying model, the Hubbard model 
with easy axis, is introduced and its validity is discussed. 
Then the dynamical mean-field theory which is used to 
investigate the correlation problem all the way from weak 
to strong coupling, as well as the Quantum Monte Carlo 
techniques employed to solve the coupled self-consistency 
equations, are discussed in Sec. [II. The results obtained 



for a half filled band and for finite doping are presented 
in Sec. IV and Sec. [v|, respectively. A discussion of the 
results in Sec. |v| closes the presentation. 
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FIG. 1. Schematic phase diagram, magnetic field H 
vs. temperature T, for a) a typical Ising-type metamagnet 
(TCP: tricritical point), b) the Ising model (|lj) in mean-field 
theory with R < 3/5 (CE: critical endpoint, BCE: bicritical 
endpoint). Solid lines: first order transition, broken lines: sec- 
ond order transition; AF: antiferromagnetic phase, P: param- 
agnetic phase. 



II. HUBBARD MODEL WITH EASY AXIS 

The Hubbard model |2(| is the generic microscopic 
model for itinerant and localized antiferromagnetism in 
correlated electron systems. For nearest neighbor hop- 
ping of electrons in the presence of a magnetic field it 
has the form 



h = -t + U J2 

NN,a i 

+ (rH)n ia , 



(2) 



where operators carry a hat. From the exact, analytic 
solution in dimensions d = 1 the (paramagnetic) ground 
state of this model at half filling (n = 1) is known to 
exhibit metamagnetic behavior, i.e. d\/dH > 0, up to 
saturation |27|. However, this is entirely due to the con- 
vex shape of the density of states in d = 1 and occurs 
even at U = 0. 

In general the Hubbard model in the form of (||) cannot 
describe metamagnetic behavior. The Hubbard model is 
isotropic in spin space and hence the antiferromagnetic 
phase described by it is isotropic, too. Consequently, any 
finite magnetic field H will orient the staggered magneti- 
zation m st perpendicular to itself. Real antiferromagnets 
are, however, never isotropic in spin space: the relativis- 
tic spin-orbit interaction L • S, where L and S are opera- 
tors for the orbital angular momentum and spin, respec- 
tively, transfers the anisotropy of position space (caused 
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by the broken rotational symmetry of the crystal lat- 
tice) into spin space, producing one, or more, easy axes 
which constrain the spins p8| . In an external magnetic 
field such a constraint leads to metamagnetic transitions, 
with or without a spin-flop depending on the strength of 
the spin-orbit interaction, as explained in Sec. |[ 

A microscopic theory of strongly anisotropic antiferro- 
magnets should ultimately be able to take into account 
the orbital degeneracy of the electrons and, by including 
the relativistic spin-orbit interaction L • S in the Hamil- 
tonian, to generate an anisotropy axis within the model 
itself. At present, this is technically not possible p9[ . 
Therefore we take the existence of the anisotropy axis e 
for granted: we employ the Hubbard model (||) and con- 
strain the magnetic moments to lie along e || H. By this 
approach the kinetic energy and the Coulomb interaction 
are treated microscopically, whereas the relativistic cor- 
rections are not. We note that the relativistic corrections 
are of the order of 1CP 2 eV and are thus small compared 
to energies of the order of 1 eV for kinetic and Coulomb 
energy. Therefore the existence of e and the correlation 
physics described by the Hubbard model (|^) are quite un- 
related. This justifies our approach where the existence 
of e is a priori assumed |30(]. 



III. DYNAMICAL MEAN- 

(D -> oo) 



FIELD THEORY 



For classical spin models (e.g. the Ising model) it is well 
known that the Weiss molecular field theory becomes ex- 
act in the limit of high spatial dimensions (d — oo). For 
lattice electrons this limit was introduced only recently 
pl| . With the proper scaling of the hopping element 
in @, t = t* j\[Z {Z = number of nearest neighbors), 
it leads to a quantum mechanical dynamical mean-field 
theory (DMFT); for reviews see |32j , |33| . The interacting 
lattice model then reduces to a self-consistent single site 
problem of electrons in an effective medium |34j] , which 
may be described by a complex, frequency dependent 
(i.e. dynamical) self energy £ CT (a->) . This problem is, in 
fact, equivalent to an Anderson impurity model comple- 
mented by a self-consistency condition |35|-|37|l. 

There are two limits in which the DMFT recovers well- 
known static mean-field theories: 

1. Weak coupling: In this situation the effective 
medium may be approximated by a static field 
which is generated by the averaged densities of the 
electrons. This leads to the Hartree-Fock approx- 
imation, e. g. E cr {ijj) — Un^ a in the homogeneous 
case, which is expected to give the qualitatively cor- 
rect behavior at weak coupling. The averaged den- 
sities n a have to be determined self-consistently. 

2. Strong coupling at half filling: Here model (^J) can 
be mapped onto the antiferromagnetic spin 1/2 



Heisenberg model for which the limit d = oo be- 
comes equivalent to the Weiss molecular field the- 
ory. 

The results obtained in these two limits will be pre- 
sented, and compared to the results for intermediate cou- 
pling, in Sec. |y[ 

The DMFT has recently provided valuable insight into 
the physics of strongly correlated electron systems, e.g. 
the Mott-Hubbard transition fl33| , f36p§| ] and transport 
properties p9| . The effect of the magnetic field H in 
(||) for n = 1 was also studied pp|-j42]|. In particular, 
Laloux et al. |^(| thoroughly investigated the magnetiza- 
tion behavior of the paramagnetic phase, assuming the 
AF order to be suppressed. For U = 3\/2i* they find a 
first order metamagnetic transition between the strongly 
correlated metal and the Mott insulator at a critical field 
H ~ 0.2t*. Giesekus and Brandt j|2| took into account 
the AF order. They considered the isotropic case where 
the field orients the staggered magnetization perpendicu- 
lar to itself, such that a metamagnetic transition cannot 
occur. 

To investigate the metamagnetic phase transition from 
an antiferromagnet to a paramagnet we consider a bi- 
partite (A-B) lattice and allow for symmetry breaking 
with respect to spin a 6 {Til} = } an d sublat- 

tice a S {A, B} = {+,-}. The self energy = 
S^(iw n ), with Matsubara frequencies u n — TtT(2n + 1), 
n = 0, ±1, ±2, . . . (using the convention K = feg = 1), and 
the Green function G a an are determined self-consistently 
by two sets of coupled equations |3l]3q] : 



G a 



de- 



N°(e) 



e 2 /z° 



(3) 
(4) 



Here z° n — iu> n + fx — and the thermal average of 
some operator 0[ip, ip*] in (§) is defined as a functional 
integral over the Grassmann variables if), ip*, with 



(Oh 



V[<il>}V[ip*}0[ili,ip*]e A <*^" '> s ' G l, (5) 



in terms of the single site action 



-U / dT^|(r)V» T (T)^(r)Vi(T) , (6) 
Jo 

where Z a is the partition function, and N°(e) is the 
density of states (DOS) of the non-interacting elec- 
trons. As the results do not much depend on its precise 
form we choose a half-elliptic DOS 2V°(e) = [(2<*) 2 - 
e 2 ]V 2 /(27ri* 2 ). cons ^ ram ^ ms( || jj j s enforced by 

setting the off-diagonal (in spin space) elements of the 
Green function equal to zero: G^} = — < VdV'I* >= 0. 
From now on t* = 1 will set our energy scale, i.e. the 
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total band width of N°(e) is equal to 4. The Dyson 
equation (||) introduces the lattice into the problem. It 
couples A— and B— sublattices and can be solved by a 
simple integration, even analytically for the above DOS. 
The functional integral (^) however is highly non-trivial 
since it couples all Matsubara frequencies. Georges and 
Kotliar |}5| and Jarrell |36| realized that the action (||) 
is equivalent to that of an Anderson impurity model, 
and can therefore be treated by standard techniques de- 
veloped for this model. Here we employ a finite tem- 
perature, auxiliary field Quantum Monte Carlo (QMC) 
method In this approach the electron-electron 

interaction is formally replaced by an interaction of in- 
dependent electrons with a dynamical, auxiliary field of 
Ising-type spins. To this end the interval [0,/3] is dis- 
cretized into A steps of size At = (3/ A. Equivalently, 
there is a high energy cut-off of Matsubara frequencies, 
i.e. = nT\2n + 1| < tt/At, n = -A/2, . . . , A/2 - 1. 
All quantities have to be extrapolated to At — > 0. The 
computer time grows like A 3 oc /3 3 , restricting A to val- 
ues below ~ 150 and (3 < 50 ... 70 on present supercom- 
puters. For small A (A < 20) one can perform a full 
enumeration (instead of the Monte Carlo sampling) of 
all 2 A possible configurations of the auxiliary field. We 
never encountered a minus-sign problem, hence no fur- 
ther approximations (like the "fix-node" method) were 
necessary. 

The self-consistency is obtained iteratively as follows: 
the Green function G (omitting indices) is calculated 
from some initial self energy, e.g. £ = 0, by the Dyson 
equation (||) . Now the new Green function G now is deter- 
mined by solving (Q) with the QMC method. Finally, the 
calculation of the new self energy E ncw = S— G~£ r +G 
completes one iteration. To improve convergence in the 
symmetry broken case G a and Gb are updated by the 
Dyson equation (^|) after every QMC simulation for one 
sublattice (Eq. [|). In the symmetry broken phases, typ- 
ically 10 — 20 iterations with 20000 MC sweeps are nec- 
essary to obtain a convergence of ~ 10~ 3 . The calcula- 
tion of a magnetization curve at (3 = 50 takes about 100 
hours on a Cray-Y-MP. Close to a phase transition the 
convergence is much slower and the statistical errors are 
larger due to strong fluctuations, in particular in the case 
of a second order phase transition. These effects limit 
the accuracy in the determination of the critical values 
of the m odel p arameters, e.g. the critical magnetic field 
(see Sec. |IVC| ). At large U- values (U > 4) the Monte 
Carlo sampling becomes more and more inefficient due 
to "sticking" problems, i.e. there are two (or more) min- 
ima in the free energy and the single spin-flip algorithm 
is no longer able to transfer between them. 

From the resulting Green functions we calculate the 
densities and the homogeneous and antiferromagnetic 
magnetization: 



(7) 



Since we are interested in the question whether the 
system is insulating or metallic we also determine the 
electronic compressibility n e = dn/dfj,. It was calculated 
both by numerical differentiation of n(/i) and from the 
two-particle correlation functions (see p4| and Appendix 
). Both results agree within the statistical errors: how- 
ever, the latter method is much more time consuming. 



IV. RESULTS FOR HALF FILLING 

A. Weak coupling 

For weak coupling ([/ < i) we expect the Hartree-Fock 
approximation to give an, at least qualitatively |l5| , cor- 
rect picture of the metamagnetic phase transitions, espe- 
cially of their order. Within this static mean-field theory 
electronic correlations are neglected and the interaction 
is decoupled as 



HF 



flifj "^fli — a ^> -r <C Titer ^> il>i—(T 

- <hi~ a xh la > . (8) 



Note, that while for e || m st || m the Fock term van- 
ishes, < cj^c^ >= 0, its presence is essential in the case 
m st _L m (see Appendix [b]). We confine our investiga- 
tions to states with homogeneous sublattice magnetiza- 
tion as described by the Ansatz 



< n 



> = 



1 



(n + a m + a a m st ). 



(9) 



Applying this Hartree-Fock decoupling scheme one ob- 
tains the effective, one-particle Hamiltonian 



N N a 

S 2 



l u 

"2~7 



— (n - am - aam st )n ia - (fj,+ <jH)h la 
2 - (m + am st ) 2 ] . (10) 



E 



This Hamiltonian is diagonalized and the one-particle 
energies e a are calculated as 



sgn(e) 




rrlZm + H). (11) 



where sgn(e) denotes the sign of the non-interacting elec- 
tron energy e. In the antiferromagnetic phase the DOS 
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has a gap of width U m st with square root singularities 
at its edge. 

From the one-particle energies ( |ll| ) the grand potential 
per lattice site (L being the number of lattice sites) is 
obtained directly as (fi, = fi — i Un) 

n/L = --InZ 

P 

= -\Hj deiV°(e)ln(l + e-^-' i )) 



U, 2 2n Un 2 
+ ^( m st+ m ) —■ 



(12) 



The potential Q has two shallow minima, one correspond- 
ing to the paramagnetic state (m ^ 0, m st — 0) and the 
other to the antiferromagnetic state (m rs 0,m st =/= 0). 
By applying a sufficiently strong magnetic field to the 
antiferromagnetic state, the paramagnetic minimum be- 
comes the lowest, such that a first order phase transition 
takes place. The reason for the occurrence of a first order 
phase transition is that, within the (static) Hartree-Fock 
approximation, the [/—term becomes minimal for largest 
(static) local moments. Therefore pure antiferromagnetic 
and ferromagnetic order are both energetically favored. 
Mixed states with m/:0 and m st 7^ - which would oc- 
cur in the case of second order phase transitions - have a 
small local moment, i.e. a high Hartree-Fock energy, on 
every second site. 

To calculate the magnetization curves the self-con- 
sistent Hartree-Fock equations are obtained from the 
minimization conditions dtt/dm st — and dfl/dm = 0: 



(13) 



™>st = — > / deiV (e — 37 



/3(e„-A) ' 



(14) 



These equations are solved numerically by iteration and 
integration according to Newton-Cotes rules. 

Within the Hartree-Fock approximation the meta- 
magnetic phase transition is found to be of first order for 
all T, even for U = 4 ( = band width); see Fig. || and [|. 
Hence a tricritical point never occurs. In this parameter 
range the Quantum Monte Carlo calculations, however, 
already show second ord er tra nsitions in a broad range 
of temperatures (see Sec. IV C ). Hence the Hartree-Fock 
solution can neither describe the experimental situation, 
where tricritical points are known to occur, nor the cor- 
rect behavior of the model for intermediate values of U. 

To estimate the anisotropy energy associated with the 
easy axis we compare the Hartree-Fock energies of the 
configurations with m st || m and m st _L m (for de- 
tails see Appendix |b]). At half filling and for U equal 
to the band width, the difference between the free en- 
ergy of these configurations does not exceed a few per- 
cent of the band width, i.e. O(10~ 2 eV). In this situation 
the spin-orbit interaction, which can be relatively strong, 



0(10 1 eV r ), indeed leads to a strong anisotropy, i.e. an 
easy axis e, along which m st is rigidly fixed. 
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FIG. 2. a) Magnetization m vs. magnetic field H in 
Hartree-Fock approximation for U — 4 at different tempera- 
tures T showing metamagnetic behavior, b) Order parameter 
for the metamagnetic phase transition (the staggered magne- 
tization) m st vs. H. The first order phase transition is clearly 



Metamagnetic phase transitions in itinerant, metallic 
systems where hitherto described by the theory of "itin- 
erant electron metamagnetism" (IEM) . In the case of an 
antiferromagnetic system in a magnetic field Moriya and 
Usami [E5| proposed a Landau theory with free energy 



F(m, m st ) 



I 



1 



2 st + am 4 + a' m\ t 



2x™ 2xst 
+b m 2 m 2 st + b' (m-m si ) 2 - Hm, (15) 

where Xm and \st are the homogeneous and staggered 
susceptibility, respectively, and the coefficients a, a', 
6, and b' are the fourth-order derivatives of the non- 
interacting free energy. Within the IEM theory the 
Coulomb interaction is treated in random phase approx- 
imation. The corresponding susceptibilities are given as 



1 



-U, 



I 

Xst 



Xst 



(16) 
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where Xm ano - Xst are ^ ne respective susceptibilities of the 
non-interacting system. The random phase approxima- 
tion for these susceptibilities is equivalent to the Hartree- 
Fock scheme described above. Therefore we may ask 
whether we obtain the IEM in the limit U <C t (where 
to, m st <C 1). The answer is not straightforward since the 
prefactors in the expansion ( |l5| ) depend, for example, on 
the lattice structure. On bipartite lattices, as discussed 
here, they diverge for T — > 0. Thus an expansion of the 
free energy in powers of m s t and m as assumed in the 
IEM Landau theory is not possible in general (for details 
see Appendix |(|). 
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FIG. 3. H — T phase diagram for different values of U in 
Hartree-Fock approximation. All phase transitions are of first 
order. Below the curves the antiferromagnetic phase is stable. 



B. Strong Coupling 

In the limit U t the Hubbard model at half filling 
(n = 1) is equivalent to an effective Heisenberg spin- 
model 



#Heis = 



J 



(17) 



NN 



where the antiferromagnetic exchange coupling is ob- 
tained in second order perturbation theory as J = 4t 2 /U. 
Spin operators are defined as Sf = gC^it — &f — 

Whkl + £ liC iT ), and Sf = -|(cj T ^ " cl ; c iT ). For this 
model the Weiss molecular field theory becomes exact in 
d = oo yielding, under the constraint of uniaxial mag- 
netization, the same results as for the Ising model. For 
Ising models metamagnetic phase transitions are well- 
studied Jul. In the case of a purely antiferromagnetic 
next-neighbor coupling (see (0)) the phase transitions 
are of first order only at T — 0, but of second order at 



all T > 0. The transition line in the H — T phase dia- 
gram has indeed the form shown in Fig. la, but with a 
tricritical temperature of T t = 0. This behavior can be 
understood already within Weiss molecular field theory, 
where the ground state energy per site is 



J* 

E(m, m st ) = -^-(to 2 - m 2 st ) - Hm, 



with 



J* = ZJ = 



IT 



(18) 



(19) 



Minimization with respect to to and m st shows that the 
fully polarized antiferromagnet (m s t = 1) has lowest en- 
ergy for H < J*/4, whereas the fully polarized ferromag- 
net (to = 1) is energetically favored for H > J*/4. Thus, 
by applying a magnetic field a first order transition is 
induced. At H = J* /A the states are highly degenerated 
since all magnetic phases with m+m s t = 1 have the same 
energy. For T > this degeneracy is lifted by entropy 
which disfavors fully polarized phases. Therefore the first 
order transition at T = immediately becomes second 
order for T > 0, i.e. T t = 0. Indeed, a tricritical point at 
a finite temperature is only obtained in the case of spin 
interactions which simultaneously favor both fully polar- 
ized antiferromagnetic and ferromagnetic configuration. 
In particular, adding a ferromagnetic interaction J' be- 
tween next nearest neighbors (NNN) on a simple cubic 
lattice as in (01) stabilizes both ferro- and antiferromag- 
netic order 

While in the case of effective spin models a ferromag- 
netic NNN coupling term is introduced ad hoc, simply 
to obtain the first order phase transition, this term natu- 
rally arises if we expand the strong coupling perturbation 
series of the Hubbard model to C(i 4 /C/ 3 ). However, be- 
sides this J'-Term there also appear additional /cmr-spin 
terms. For the hypercubic lattice the effective Hamilto- 
nian H c g reads pq| 



NN 

t A 

—A (Q12Q34 + Q14Q23 - Q13Q24) , (20) 



u u + B W 

t 4 



J ' = a Tp c ' 



(21) 



(22) 



Here r and r' are lattice vectors connecting a site to its Z 
neighbors, and □ represents a plaquette. Each plaquette 
is counted only once: the four sites {1,2,3,4} represent 
its four corners in clockwise or anticlockwise order. The 
constants A, B and C depend on the lattice, and the 
Hermitian operators Qij are defined as 



Qi 



2 s, ■ s„ 



(23) 



G 



The plaquette contribution competes with the ferro- 
magnetic NNN term (J') and drives the system to sec- 
ond order phase transitions. For the hypercubic lattice 
the plaquette contribution is stronger than the ferromag- 
netic NNN term, yielding second order phase transitions 
even for T = 0. The same is true for the Bethe lattice 
where, in fact, J' < (for details see Appendix ^). Thus 
in strong coupling perturbation theory the metamagnetic 
phase transition is of second order even at T = 0. 

C. Intermediate coupling 

The perturbation analysis described above demon- 
strates that the order of the metamagnetic phase tran- 
sition depends on the Coulomb interaction U in a deli- 
cate way. For small U the phase transition is purely of 
first order and for large U of second order. Apparently, 
the tri- or multicritical point linking these two regimes 
must be found at intermediate coupling. In this im- 
portant, non-perturbative regime Quantum Monte Carlo 
techniques are employed to solve the problem numeri- 
cally without any further approximation. The results 
for the magnetization m(H) and the staggered magne- 
tization m s t(H) are shown in Fig. |] for U = 2. Below 
the Neel temperature a metamagnetic behavior is clearly 
seen: for small magnetic fields the magnetization is ex- 
ponentially suppressed with temperature. Then, towards 
the metamagnetic phase transition, the susceptibility in- 
creases drastically and becomes maximal at the critical 
field H c . Second order phase transitions are observed for 
1/14 <T <T N = 0.114 ± 0.006, whereas the transition 
is of first order at lower temperatures, i.e. T < 1/16. At 
the phase transition the order parameter, i.e. the stag- 
gered magnetization, vanishes. From the curve m s t(H) 
the critical field H c and also the order of the phase tran- 
sition is determined by a square root fit for second order 
transitions and by the mean of the hysteresis for first 
order transitions. 

Using these values of H c the phase diagrams, Fig. [jj, 
for different values of U are constructed. The case U = 4 
( = band width) and half filling, Fig. |^a, was already 
presented in Ref. fi^ . This phase diagram shows both 
first order (for T < 1/16) and second order phase transi- 
tions (for 1/8 < T < T N « 0.2). The field dependence at 
intermediate temperatures, i.e. 1/16 < T < 1/8, is more 
complex: m st is almost field-independent in region AFi 
(where m ~ 0), decreases sharply at the boundary to re- 
gion AFn (where m > 0) and, upon further increase of 
the field, vanishes in a second order transition. Although 
the error bars do not permit an unambiguous interpre- 
tation it hence seems that the order parameter decreases 
by two consecutive transitions: the first one being of first 
order or corresponding to an anomaly, and the second 
one being of second order. Taken together the results 
seem to correspond to the scenario of Fig. @b. For U = 2 
phase transitions of first order are found for T < 1/16 



and of second order for 1/14 < T < T N = 0.114 ± 0.006 
(see Fig. ||b). Here, the temperature regime with two 
consecutive transitions, as obtained for U = 4, has disap- 
peared or has become very small: the scenario is similar 
to Fig. Ilia. The phase diagram for U — 3, also displayed 
in Fig. |^b, shows some features of the phase diagram in 
Fig. |l|b or U = 4, respectively. In particular one observes 
a maximum in the second order phase transition line and 
different slopes for the second and first order line at the 
crossover point. However, the numerical data do not in- 
dicate the existence of two consecutive transitions. Thus 
this phase diagram lies in between the scenarios depicted 
in Figs. |l|a and [l]b. 




H 

FIG. 4. QMC results, including error bars, a) for the mag- 
netization m(H) as obtained for the d = oo Hubbard model 
with easy axis at half filling and U = 2 (= half band width) for 
different temperatures, b) Staggered magnetization mst(H) 
for the two temperatures below Tm- 
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FIG. 5. H — T phase diagram for the d — oo Hubbard 
model with easy axis at half filling as constructed from the 
QMC results for m(H) and m at (H), a) U = 4, b) U = 2,3. 
Second order phase transitions are indicated by dashed lines, 
first order transitions by solid lines. Curves are guides to the 
eye only. 



To study the influence of the interaction U qualita- 
tively and quantitatively, it would be desirable to calcu- 
late H — T phase diagrams at even larger values of U . 
Unfortunately the Quantum Monte Carlo approach fails 
in this c ase due to the problem of "sticking" as mentioned 
in Sec. 



[II 



Therefore, we concentrate on the crossover 
from the intermediate coupling regime with first and sec- 
ond order transitions to the weak coupling regime with 
first order transitions only. The results for half filling 
are collected in Fig. ^ showing the [/-dependence of two 
transition temperatures: the top curve is the Neel tem- 
perature T N at H = (taken from Ref. Q). The lower 
curve corresponds to the temperature T c where the sec- 
ond order phase transition line terminates, i.e. it rep- 
resents either the tricritical or the critical temperature 
of Fig.|l]. For temperatures below T c a first order meta- 
magnetic phase transition is observed in an external mag- 
netic field. Fig. H reveals the crossover from intermediate 
coupling with first and second order phase transitions to 



weak coupling with first order transitions only: as U de- 
creases the regime with second order phase transitions 
(T c < T < Tjy) shrinks, while the temperature regime 
for first order transitions remains nearly unchanged up 
to U re 2. 
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FIG. 6. Neel temperature Tjv (circles) and (tri-)critical 
temperature T c (diamonds) vs. U. Above Tjv the system is 
paramagnetic. In an external magnetic field the order param- 
eter vanishes in a second order metamagnetic phase transition 
for T c < T < Tn and in a first order transition for T < T c , 
respectively. 



V. RESULTS AWAY FROM HALF FILLING 

In the preceding sections metamagnetic transitions 
were investigated in the case of half filling. Beyond 
half filling the commensurate antiferromagnetic phase re- 
mains stable in the parameter regime under considera- 
tion ( 6 = |1 - n\ < 0.075, T > 1/32). Incommensurate 
spin density waves become stable only in a small density 
regime at lower temperatures f4§|| . Another possible in- 
stability of the antiferromagnetic phase away from half 
filling is phase separation, found within the Hartree-Fock 
approximation and in second order perturbation theory 
at constant order parameter at T=0 [^9|. However, at 
least for T > 1/16 we do not observe phase separation 
since the electronic compressibility n e = dn/dfi is finite 
and positive (see Fig. 0). 

Upon doping the magnetization curve changes consid- 
erably and hardly indicates the existence of a metamag- 
netic phase transition (Fig. ||). This is due to the fact 
that in the metallic phase there is no longer a "Slater 
gap" at the Fermi energy; therefore the homogeneous 
susceptibility is not as strongly affected by the antiferro- 
magnetic order as at half filling. The phase transition is, 
however, clearly seen in m st (H). 
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FIG. 7. Particle number n vs. chemical potential fi as cal- 
culated by grand canonical QMC simulations in the antifer- 
romagnetic phase for U = 4. 
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FIG. 8. Change of a) the magnetization m(H), and b) the 
staggered magnetization m s t(H) with doping for U = 2 at 
T = 1/16. 



From the m s t vs. H curve the phase diagram (Fig. g) 
is constructed. The metamagnetic phase transition line 
is found at lower temperatures and fields compared to 
half filling. 

Associated with the metamagnetic phase transition is 
a change of the electrical resistivity. To study this im- 
portant effect we calculated the ^/-dependence of the 
electronic compressibility K e . This quantity indicates 
whether the system is metallic or insulating. For an insu- 
lator n e vanishes for T = and is exponentially small for 
temperatures lower than the antiferromagnetic gap. In a 
Fermi liquid, on the other hand, n e (T — 0) is finite since 
it is proportional to the density of states at the Fermi 
level and hence proportional to the Drude conductivity. 
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FIG. 9. Change of the H — T phase diagram with dop- 
ing 8 for [7 = 4. QMC data are shown for a) 8 = 0, b) 
8 = 0.025, 0.05, and 0.075. Second order phase transitions 
are indicated by dashed lines, first order transitions by solid 
lines. For T = 1/32 the numerical error does not permit the 
determination of the order of the phase transition unambigu- 
ously. 



9 



The results for n e as a function of magnetic field H at 
U = 2 are shown in Fig. [H| At half filling, 5 = 0, the 
compressibility is seen to increase with H . This effect is 
particular pronounced at low temperatures (T = 1/25) 
where K e is essentially zero at low fields and rises to 
n e w 0.3 above the critical field, indicated by an ar- 
row. Hence the metamagnetic phase transition is a tran- 
sition from an antiferromagnetic insulator to a metal 
with homogeneous magnetization. At higher tempera- 
tures, T = 1/14, the compressibility is always finite due 
to thermal excitations. We note that at U = 4, when 
the electrons are essentially localized, n e remains small 
(0 < n e < 0.03 at T = 1/8; not shown in Fig. |HJ) 
even above the critical field, indicating an insulator-to- 
insulator transition. The situation is very different at 
finite doping (S = 0.05). Here Fig. 10 shows that K e de- 
creases with H by approximately 50% as the system goes 
through the metamagnetic transition from an antiferro- 
magnetic metal to a metal with homogeneous magneti- 
zation. 
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FIG. 10. Field dependence of the electronic compressibility 
K e = dn/d/j, for U = 2 at S = (T = 1/14, 1/25) and S = 0.05 
(T = 1/14). The arrows mark the respective critical fields for 
the metamagnetic phase transition. 



VI. DISCUSSION 

In summary, we investigated the origin of metamag- 
netism in strongly anisotropic antiferromagnets starting 
from a microscopic model of strongly correlated elec- 
trons, the Hubbard model with easy axis, by employing 
a dynamical mean-field theory. This approach is fun- 
damentally different from previous investigations since 
we identified and explicitly evaluated the simplest elec- 
tronic, i.e. fully quantum mechanical, correlation model 
that is able to explain the conditions for metamagnetism. 
For this electronic model we show unambiguously that at 
intermediate coupling the phase transition is of first or- 



der at low temperatures and of second order near the 
Neel temperature, i.e. the order of the phase transition 
changes. This phenomenon, which has been of interest to 
various communities in classical statistical mechanics for 
a long time already, is extracted from a model of itinerant 
electrons. 

Our approach allows us to describe a broad range 
of qualitatively different metamagnets within a single 
model. While at present this simple model does not per- 
mit any quantitative calculation of material properties 
it does describe itinerant and localized, metallic and in- 
sulating metamagnets and the crossover between them. 
This crossover is related to two fundamental experimen- 
tal parameters, i.e. pressure (related to U/t which de- 
creases with pressure) and doping. 

At half filling the Coulomb interaction leads to a 
crossover from a band insulator to an insulator with lo- 
calized moments. Thereby the phase transition changes 
from first order for the bandlike metamagnet to second 
order for the localized one. Only at intermediate cou- 
plings are both first and second order phase transitions 
observed as found in experiment. The H — T phase dia- 
gram obtained for an intermediate Coulomb interaction 
{U = 4 = band width) is strikingly similar to that of 
FeBr 2 @f(| or the Ising model with J' < J @,|l9|. 
We note that in these insulating systems the applicabil- 
ity of a theory which becomes exact for large coordination 
number is justified by the fact that the AF superexchangc 
involves 20 equivalent sites in the two neighboring iron 
planes |S1[ . At smaller values of the Coulomb interaction 
(U = 2) the temperature regime with second order transi- 
tions shrinks and the two step phase transition becomes 
less pronounced, reproducing the scenario of Fig. la, as 
observed e.g in FeCl2 ||. 

The calculations off half filling allow us to investi- 
gate the properties of metallic metamagnets, such as the 
Uranium-based mixed-systems |^|, for which a theory 
in terms of a correlated electron model is mandatory. In 
contrast to the insulating case, the metamagnetic phase 
transition in the metallic system is hardly visible in the 
magnetization curve. This is because there is no longer 
a gap at the Fermi energy. Quite generally, the critical 
temperatures and fields decrease upon doping. 

The metamagnetic transition is accompanied by pro- 
nounced changes in the conductivity of the system. The 
Hubbard model with easy axis can qualitatively describe 
several scenarios: 

(i) In the insulating, localized regime (U > 4 at half 
filling) a magnetic field causes a transition from an an- 
tiferromagnetic insulator to an insulator with homoge- 
neous magnetization. 

(ii) At lower [/-values (e.g. U = 2) at half fill- 
ing an insulator-to-metal transition occurs at the mag- 
netic field where the AF order disappears. Such a phe- 
nomenon is observed, for example, in the AF phase 
of Lai_ x Caa;Mn03, where the resistivity is found to 
change by several orders of magnitude J5*|. This is re- 
ferred to as "colossal" magnetoresistance. We note that 
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Lai_a;Ca x Mn03 shows no strong anisotropy. Therefore 
our approach can only describe the general features, in 
particular the existence of the insulator-to-metal transi- 
tion. 

(iii) Away from half filling a magnetic field induces 
a transition from a metallic antiferromagnet to a metal 
without staggered moment. Here the compressibility 
changes by less than an order of magnitude, e.g. about 
50% at U = 2, 6 = 0.05. A similar effect is found in sev- 
eral strongly anisotropic antiferromagnets, both in multi- 
layers and bulk intermetallic compounds such as UPdGe 
H . In these systems the origin of this "giant" magne- 
toresistance is attributed to band structure effects and 
spin scattering |53|]. By contrast, our approach stresses 
the importance of genuine electronic correlation effects. 

More detailed investigations, including band degener- 
acy and spin-orbit interaction, may eventually provide 
even quantitative insight into these interesting and im- 
portant phenomena. 



Acknowledgments 

We acknowledge useful correspondence with N. Gior- 
dano, G. Lander, B. Liithi, A. Ramirez, and W. Wolf, 
and are grateful to H. Capellmann, P. van Dongen, W. 
Mctzner, H. Muller-Krumbhaar, J. Schlipf, F. Steglich, 
G. Stewart and, in particular, to V. Dohm, W. Klecmann 
and W. Selke for very helpful discussions. In its initial 
stages this work was supported in part by the SFB 341 
of the DFG and by a grant from the ONR, N000 14-93- 
1-0495. 



dx 



rp \ a a 

/ j nn' ,n' 



n l an ' 



(A3) 



where T is the local two particle correlation function 
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The quantity t£* = ^{{G^ 1 + ££„} in © mea- 
sures the response of the averaged medium to an infinites- 
imal change of the field x. This dynamical response func- 
tion is determined by an integral equation in frequency 
space which does not explicitly depend on momentum. 
(Note that there are no convolutions in k space in the 
d = oo limit as is typical for a mean-field theory). This 
property does not imply, however, that the response func- 
tion 7^ is local, too. It only indicates that 7^ is diag- 
onal in the momentum k. Momentum dependence enters 
implicitly by the particular k dependence of the exter- 
nal field (k — in the case of the compressibility or the 
ferromagnetic susceptibility, and k = (tt, . . . , it) for the 
staggered susceptibility). 

In the presence of an external field the variables z° n 
in the Dyson equation (|J) are replaced by z° n = iuj n + 
fi + <jH + aaH st + a^cDW — £w- The derivative of the 
Green function yields: 
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APPENDIX A: CALCULATION OF 
SUSCEPTIBILITIES FROM CORRELATION 
FUNCTIONS 

Quite generally susceptibilities can be obtained from 
the derivative of the order parameter m x w.r.t. the cor- 
responding field x: 



Xx 



dm x 
dx 
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deD(e)(z° an -ey Z °_ an )- 2 , (A6) 
-OO 

/oo 
deD(e)(z: n -ey Z ° an )- 2 e 2 /{z°_ an ) 2 . (A7) 
-oo 

Since ( |A5| ) separates in Matsubara frequencies n and spin 
a it can be easily solved for dG/dx: 



dx 



EB» aa I ra _ ox \ 
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with r ax = 




for x = 




(A2) 



Here x = H and x — H s t lead to the ferromagnetic and 
antiferromagnetic susceptibilities, x = [i to the electronic 
compressibility, and x = ficDW to the charge density 
wave susceptibility. From the two self-consistency equa- 
tions (||) and (|]) one obtains two corresponding equations 
for the derivative of the Green function w.r.t. the variable 
x. The derivative of the functional integral (H) gives: 



with the 2x2 array in a G {^4, B}: 

RS = -(detD)- 1 D^ ) 
whereby D and T are defined as 
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Now dG/dx can be eliminated by setting ( |A3| ) equal 
to (|A^), yielding 

Ja'x 

a' 

]T {M^^^+^TC^J^J. (A12) 

From this equation we determine 7 by numerical inver- 
sion of a 4A x 4 A matrix. Knowing 7 we obt ain dG/dx 
via (A3) or and thus the susceptibility (Al). 



APPENDIX B: HARTREE-FOCK THEORY FOR 

m st -L m 

Similar to the derivation of the Hartree-Fock equations 
for m || H || m st (sec Sec. IVA) we will now investigate 
the case with perpendicular orientation m || H _L m st . 
The Ansatz for the one particle densities 



1 



1 , N 

s >= -am st (Bl) 



<n ieaa >= -{n + am),<c iea ^ iE(x9 . } 

yields in addition to the Hartree term a Fock term in the 
decoupling (p|) 
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With this Ansatz one readily obtains the effective one- 
particle Hamiltonian 

w HF = J2 tiA^ - \\ E(" 2 - m2 - m2 *) + ( B3 ) 
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-(n - am)n ia - —am s tC[,f m 



Diagonalizing this Hamiltonian yields the one-particle 
energies 
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From these energies the free energy ft is calculated, and 
the minimization with respect to m and m s t leads to the 
following Hartree-Fock self-consistency equations 
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APPENDIX C: SERIES EXPANSION OF THE 
HARTREE-FOCK FREE ENERGY 

The itinerant electron metamagnetism theory of Mori- 
ya and Usami |2f| can be derived from the Hartree-Fock 
approximation only if the free energy is analytic in m 
and m s t- Since in RPA the Hubbard interaction U con- 
tributes to the free energy analytically (see ( |l5|) and (|l^)) 
any non-analytic behavior must be due to the kinetic en- 
ergy. Its expansion in the order parameter m st at T = 
is analyzed in this section. 

To calculate the expansion in m st a staggered magnetic 
field H s t is introduced: 



NN,a 
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On A-B lattices the one-particle energies for this Hamil- 
tonian show a gap at e = with square root singularities 
at its edge 



e = sgn(e) x /e 2 + J ff s 2 t . 



(C2) 



We consider the half filled band, where the staggered 
magnetization m st is calculated from the one-particle en- 
ergies. It shows the following asymptotic non-analytic 
behavior for H st — > 0: 



-H st 



2N°(0)H st ln(l/H st ) + 0(H s 



(C3) 
(C4) 



Similarly the asymptotic behavior of the energy (CI) is 
obtained as 



AE(m s t)=J2 j deA°(e)(e~-e) 



-H st m st + 0(H 2 st ). 



(C5) 
(C6) 



Subtracting the contribution due to H st , AEn st — 
—H st m st , the asymptotic dependence of the kinetic en- 
ergy on m st reads 



AE kin (m st ) 



H 3t ^0 1 



2 ln(l/m at ) 



(C7) 



This shows that AE^i n (m st ) is non-analytic in m st . 
Therefore the itinerant electron metamagnetism theory 
cannot be derived from the Hartree-Fock theory for the 
Hubbard model with an easy axis. 
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APPENDIX D: THE METAMAGNETIC PHASE 
TRANSITION AT STRONG COUPLING 



As in Sec. [V A these Hartree-Fock equations are solved 
numerically. 



In the limit of strong coupling and half filling the 
0(t A /U 3 ) perturbation theory yields the effective spin 
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Hamiltonian (|20|). In the following, we study the meta- 
magnetic phase transition, and especially the order of 
the transition, for this effective Hamiltonian. Restrict- 
ing ourselves to solutions with mixtures of ferromagnetic 
(m) and antiferromagnetic (m s t) order, the ground state 
energy is a polynomial in m and m st (t* = 1): 



1 1 



1 1 



E = m mz f 

2U 2 U st 
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-Am 4 - 2 Am 2 m 2 st - 2 A m 2 + A m\ t + 6 A m 2 s 



-16Cm 2 - 16Cm 2 t } - i? m + const. 



(Dl) 



One can see that the ferromagnetic next nearest neighbor 
interaction C favors both saturated antiferromagnetism 
and saturated ferromagnetism rather than ferrimagnetic 
phases. By contrast the plaquette term A has contribu- 
tions that support the formation of a ferrimagnetic state. 
To obtain the ground state the energy must be minimized 
with respect to m and m st under the constraints |m| < 1 
and |m s t| < 1 — |m|. Differentiation of E with respect 
to m s t shows that, for fixed m, E has one maximum at 
m st = and two minima at 



± 



8U 2 + Am 2 -3A + 8B + 8C 



A 



For sufficiently strong coupling U (e.g. U > 2 in the case 
of the hypercubic lattice with A=20, C=2 and B=A) 
these minima are outside the constraint \m s t \ < 1 — \m\. 
Therefore E becomes minimal at the border of the con- 
straint, i.e. for \m st \ = 1 — \m\. Replacing m st by 1 — m 
the minimization with respect to m readily yields for the 
ground state 







for 



tt s 1 -A+2C+2U 2 +2B 

u ^ 2 m 



for H > - c+ $+ B 



A-2C-2U 2 -2 B+2 H IT 
A-4C 



else 



m st 



1 



This ground state solution for the effective spin Hamil- 
tonian shows a second order metamagnetic phase tran- 
sition for A > AC . This is the case for the hypercubic 
lattice (A = 20, C = 2) and for the Bethe lattice, where 
^4=0 but C = — 1, i.e. the next nearest neighbor coupling 
is antiferromagnetic. In conclusion, the strong coupling 
theory shows second order phase transition for all tem- 
peratures, even at T = 0. 
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